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We first review tiie general frameworlt wliicli enables one to resum high-energy logarithms in 
hadronic processes, both in the evolution of parton densities and in the coefficient functions. We 
then present an all-order calculation in perturbative QCD of the inclusive Drell Yan and vector boson 
production in hadron-hadron collisions, in the limit where centre of mass energy is much bigger 
than the invariant mass of the final state. Our calculation resums leading non-trivial logarithms 
in the ratio of these two scales. We also study the phenomenological relevance of the high-energy 
corrections for Drell-Yan processes at the LHC. We find that the resummation corrects the NNLO 
result by a few percent, for values of the invariant mass of the lepton pair below 100 GeV. Corrections 
to W ,Z production are expected to be of the same order. 



I. QCD IN THE HIGH-ENERGY LIMIT 

Perturbative calculations in QCD with initial-state 
hadrons are possible because of asymptotic freedom and 
coUinear factorization. In particular the latter enables us 
to separate, in the presence of a hard scale, pertburba- 
tive coefficient functions from universal parton densities. 
However, when the hard scale Q is well below the centre- 
of-mass energy of the colliding hadrons v^, the typical 
values of momentum fraction x of the colliding partons 
may be rather small: X1X2 — / S 1. Whenever x 
is small, logarithms of x may spoil the perturbative se- 
ries. Thus accurate calculations require the computation 
of the coefficients of these logarithms, and if they are 
large it may be necessary to resum them. High-energy 
logarithms contaminate both the perturbative coefficient 
functions and the evolution kernels for the parton densi- 
ties. 

The resummation of small- a; logarithms in the pertur- 
bative evolution of parton distribution functions at the 
next-to leading logarithmic accuracy (NLLz) has been 
performed by different groups (see for example Ref. [Ij 
for a comparative review). Any reliable small- a; resum- 
mation should address the issue of the perturbative in- 
stability of the BFKL kernel; the resummed anomalous 
dimension should be computed at the NLLa; accuracy 
and must match standard DGLAP at moderate values 
of X. We shall briefly review the resummation procedure 
proposed by Altarelli, Ball and Forte (ABF), which, in 
particular, provides a consistent treatment of the factor- 
ization scheme dependence. This is an essential ingredi- 
ent if one wants to match the resummation to fixed order 
calculations. For a more detailed description of the ABF 
resummation see for instance [2] and references therein. 

The general procedure for resumming inclusive hard 
cross-sections at the leading non-trivial order through 
fcT-factorization is known and its implementation 
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when the coupling runs understood j5|. High-energy re- 
summed coefficient functions are now known for an in- 
creasing number of processes. Calculations have been 
performed for photoproduction processes [H, fsj , deep in- 
elastic processes 0, |j] , hadroproduction of heavy quarks 
0) &0| ) and gluonic Higgs production both in the point- 
like limit fsl], and for finite top mass mt In par- 
ticular the latter result has been used, together with the 
l/mt expansion, to study finite top mass effects in the 
NNLO inclusive cross-section |ll| - [l3| . The resummation 
for Drell-Yan and vector boson production has been per- 
formed in |l4l| and, more recently, for direct photon pro- 
duction in [l5|. Thanks to this effort it is now possible 
to study the impact of small- a; resummation at the LHC. 

II. DRELL-YAN AND VECTOR BOSON 
PRODUCTION CROSS-SECTION AT THE LHC 

Accurate perturbative calculations of benchmark pro- 
cesses such as Drell-Yan production of a lepton pair and 
the closely related vector boson production cross-section 
are an essential component of the LHC discovery pro- 
gramme. Outstanding precision has been reached in the 
calculation of radiative corrections in perturbative QCD. 
The inclusive cross-sections have been known at NNLO 
for a long time The calculation of the rapidity 

distribution at the same accuracy has been performed 
in [2^. More recently the fully differential cross section 
has also been computed. The resummation of thresh- 
old logarithms is known up to N^LL [2ll - [2^ . Thanks 
to these very precise calculations, the theoretical uncer- 
tainty for W^/Z production is estimated to be below 4%. 
The dominant contribution comes from parton distribu- 
tion functions (1.5 — 2.5%) [25l - [27j . while higher order 
corrections are estimated by scale variations to be below 
1%. Other sources of uncertainty, such as the treatment 
of heavy flavours and small- a:; effects, are more difficult 
to assess. An estimation of the latter ones is the target 
of this study. It is important to control the theoreti- 
cal calculation at the same precision as the experimental 
accuracy, which is claimed to be around 4%, see for in- 
stance [231 , so that these processes can be used to monitor 
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the parton luminosity at the LHC. 

Moreover the LHCb experiment aims to measure the 
Dreh-Yan production of a muon pair at invariant mass 
as low as 2.5 GeV [S^l- Such a measurement will probe 
values of x down to 10~^ and will reduce the uncertainty 
of the largely unconstrained parton distribution functions 
at very small values of x. On the theoretical side the 
above mentioned fixed order calculation of the Drell-Yan 
rapidity distribution shows that the perturbative series 
is unstable for invariant mass of the muon pair below 
20 GeV, suggesting important contributions from small- 
X evolution [29|. 

III. HIGH-ENERGY RESUMMATION 




(a) (b) 

FIG. 1: Feynman diagrams for the process g*{k) qijpi) — ^ 
liPi)- The quarks are on-shell and massless, while the 
initial state gluon and the final state photon are off-shell. 



A. Resummation of the evolution kernels 

The resummation of collinear and high-energy loga- 
rithms can be achieved by considering the DGLAP and 
BFKL evolution equations respectively. The DGLAP 
equations describe the dependence of the parton den- 
sities 



^{N,as{t))G{N,t), 



(1) 



where t = In ^ and the N is the Mellin moment with 
respect to x. The high-energy, or small-x, limit then cor- 
responds to N 0. We consider the largest eigenvalue 7 
of the anomalous dimension matrix in the singlet sector, 
which has been computed to NNLO accuracy 3l|. The 
BFKL equation can be written as an evolution equation 
for the Mellin moments of the unintegrated gluon distri- 
bution 

^g{^,M) = x{M,as{dM))g{i,M) , (2) 

where ^ = In ^ and M is the Mellin moment with respect 
to Q^; therefore the collinear limit in Mellin space is the 
M ^ limit. The differential operator as{dM) repre- 
sents the running coupling in Mellin space; at one-loop 
we have 
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dM 



The BFKL kernel x has been computed explicitly at the 
next-to- leading order accuracy j32l |. The fixed-order ex- 
pansion of the kernel is known to be slowly convergent: 
the NLO order term xi is large and has a qualitatively 
different shape to the leading order kernel xo- The ori- 
gin of this instability is the behaviour of the kernel in the 
collinear limit: Xo ^ jji while xi ~ The first step 

towards a reliable evolution at small-x is the resumma- 
tion of the collinear singularities of the BFKL kernel. In 
the ABF approach this is achieved by exploiting duality 
relations among the evolution kernels: 



M = j(x(M,as),as), 
N = xh{N,as),as). 



These relations state that at fixed coupling the two evo- 
lution kernels are the inverse of one another. Therefore x 
determines the high-energy (TV = 0) singularities of the 
DGLAP anomalous dimension and 7 the collinear poles 
(Af = 0) of the BFKL kernel to all orders. Thus we can 
construct a double leading expansion of the BFKL kernel 
to next-to leading order: 



= (xo + Xs- Xo") + (Xi + Xs 



xf') 



(4) 



(3) 



where Xs and Xss are dual to the LO and NLO DGLAP 
anomalous dimension respectively; therefore collinear 
logarithms are resummed to NLL. The terms Xd c avoid 
double counting in the two expansions. The double lead- 
ing kernel, once symmetrized, is stable in the whole re- 
gion < M < 1 and can be used to compute an anoma- 
lous dimension which resums both logarithms of and 
X. However the resummation aims NLLx accuracy and 
other effects such as running coupling and scheme de- 
pendence must be addressed (see Ref [2] and references 
therein). The formalism to consistently resum small- a; 
logarithms in the coUinearly factorized coefficient func- 
tion is discussed in the next sub-section. 



B. Resummation of the coefficient function 

The leading non-trivial high-energy singularities for 
inclusive cross-sections can be obtained using the kx- 
factorization theorem [1, [I]- In [3 we have ex- 
tended such an approach to deal with the rather com- 
plicated fiavour structure which characterizes Drell-Yan 
processes. The main result of that computation is the 
resummed coefficient function for the quark-gluon sub- 
process. It can be computed by considering the process 

9*{k)q{pi)^j*{q)q{p2), 

in the framework of fcT-factorization. The relevant Feyn- 
man diagrams are shown in Fig. [U where = P2 = 0, 
but ^ -|k|2 and <f- = Q'^ . The double N,M Mellin 
transform of the off-shell cross-section, in the small TV 
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limit is found to be 



resummed kernel. We define 



a. AT{l-MfT{l + Mf 



(5) 



M)(2-Af)(3- Af) ■ 

An implicit relation for the resummed coefficient function 
Dqg can be found in terms of the above function hqg^ the 
small- a; resummed anomalous dimensions ^gg and ^qg and 
the scheme dependent factor R (Jj: 



lqg+lggDqg{N,a,) = /l,g (0, 7gg)i?(7gg) 



(6) 



The analytic expression for the leading high-energy be- 
haviour of the coefficient function can be computed to 
any desired power of the strong coupling. In Mellin space 
the leading iV — behaviour of the MS coefficient func- 
tion for the qg channel is 



Dqg{N, as) 
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The coefficients of O (a^) and O (a^) are in agreement 
with the high e nerg y limit of the fixed order NLO [l^ 
and NNLO j3 computations. This is a very non- 
trivial check of the procedure. The O (a^) and subse- 
quent terms are all new results. The high energy singu- 
larities of the quark-quark coefficient function are now 
easily deduced using the colour-charge relation: 



Dqq{N, as) 



Cf 
Ca 



Dqg{N,as)-j^TR 



(8) 



which is valid in the high energy limit. It can be shown 
that all the other coefficient functions are sub-leading in 
the high energy limit [13] • At the logarithmic accuracy 
considered in this study, the high-energy singularities of 
the vector boson production cross-sections are the same 
as those for Drell-Yan. 



IV. LHC PHENOMENOLOGY 

We now to study the impact of high-energy resumma- 
tion on the inclusive Drell-Yan cross-section at the LHC. 
We start by considering the resummed result matched 
to the fixed order calculation at 0{as) using NLO par- 
ton distribution functions from the NNPDF collabora- 
tion [3^, with the evolution performed using the ABF 
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Df^O ^ JNLO ^ JNLO 



(9) 



where Dj^'^^ are the coefficient functions for the different 

''J 

hard sub-processes and fi's are the parton distribution 
functions. This JsT-factor is plotted in Fig.[2]as a function 
of Q (top line in blue) at the LHC centre of mass energy 
VS = 14 TeV. At large Q it approaches one, as it should; 
at Q = 100 GeV the resummation corrects the fixed order 
result by 7 % and the effect is as big as 15 — 20 % at Q 
10 GeV. We then study the effect of the resummation on 
the NNLO calculation. Because the ABF resummation 
has not been implemented at this order, we consider two 
different ratios: 
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This definition of the if-factors is such that they again 
approach one at large Q; they are plotted in Fig. [21 
j^NNLO green in the middle and K^^^'^ in red, bot- 
tom curve. The effect of the resummation is diminished. 
In particular K^'^^'~' is always very close to one, sug- 
gesting the the NNLO coefficient functions capture most 
of the small- a; behaviour of the hard process. In K^^^'^ 
instead the high energy contributions to the parton evo- 
lution is also taken into account. In this case the re- 
summation corrects the fixed order results by 5 — 7% for 
Q < 100 GeV. Thanks to these two results, we can es- 
timate that the resummation corrects the NNLO by no 
more than a few percent for Q < 100 GeV. Since the 
effects are rather small, and the choice of parton distri- 
bution functions make a sizeable effect on the impact of 
the resummation, to obtain a reliable calculation of the 
resummation corrections it will be necessary to refit the 
PDFs using resummed evolution and cross-sections. This 
is in principle feasible, but it has not yet been done. 

The plot also shows that small- a; effects are potentially 
significant for the W/Z cross-section {Q ^ mz/w), where 
the theoretical uncertainty is claimed to be of the order 
3 — 4%. However more detailed studies must be per- 
formed: the LLx resummation is the same for Drell-Yan 
processes and W/Z production but the fixed order results 
do change. 

In Fig. [3] we compare the relative size of some of the 
various contributions beyond NLO. To this purpose we 
define four X-factors with the same denominator: the 
NLO cross-section computed with NLO partons. The 
two blue curves at the top are the NLO resummed cross- 
section with resummed partons (solid line, same as in 
Fig. [21) and the NLO resummed cross-section with stan- 
dard NLO partons (dashed line). The curves at the bot- 
tom are the NNLO resummed cross-section computed 
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Q/GeV 

FIG. 2: In this plot from top to bottom we have: the NLO Jf-factor as defined in Eq. ([9| and the two NNLO X-factors of 
Eq. m, and J^f '^■^o respectively. 




Q/GeV 

FIG. 3: In this plot we have the NLO (top curves) and NNLO (bottom curves) resummed cross-section, normalised to NLO. 



with NNLO parton distributions (solid green) and NLO 
ones (dashed green). The gluon-gluon subprocess makes 
a substantial contribution at NNLO. However this contri- 
bution has not been resummed yet: to do this would re- 
quire the calculation of the off-shell process g*g* 'y*qq- 



V. CONCLUSIONS 

We have presented analytic results for the Drell-Yan 
(and thus also vector boson production) coefficient func- 
tions to arbitrarily high orders in the strong coupling, 
in the limit of high partonic centre-of-mass energy. Our 
results are given in MS scheme and they agree with the 



known results at NLO and NNLO, while providing new 
results at N^LO and beyond. 

We have evaluated the effect of high-energy resum- 
mation by computing K-factors between resummed and 
standard fixed order results. We have considered proton- 
proton collision at — 14 TeV. We have found that 
the resummation corrects the NLO result by 5-10 % for 
Q < 100 GeV. When NNLO corrections are included the 
effect is reduced to a few percent in the same kinemati- 
cal region. More studies have to be done, especially for 
vector boson production, where the effects of small-x re- 
summation could potentially be of the same order of the 
current theoretical uncertainty. 
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